Spin Reduction Transition in Spin-| Random Heisenberg Chains 
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Random spin-| antiferromagnetic Heisenberg chains are investigated using an asymptotically exact 
renormalization group. Randomness is found to induce a quantum phase transition between two 
random-singlet phases. In the strong randomness phase the effective spins at low energies are Seff = |, 
while in the weak randomness phase the effective spins are Seff = |. Separating them is a quantum 
critical point near which there is a non-trivial mixture of spin-i, spin-1, and spin-| effective spins at 



low temperatures. 



Some of the most dramatic effects of randomness in 
solids appear in the low temperature behavior of quan- 
tum systems. A (deceptively) simple class of such systems 
are random quantum spin chains, in particular Heisenberg 
antiferromagnetic chains with Hamiltonian 



H — ^ JiSi ■ Si^ 



(1) 



From a real-space renormalization-group (RG) analysis, 
it has been shown that the spin-i random anti- 
ferromagnetic (AFM) chain is strongly dominated by ran- 
domness at low temperatures even when the disorder is 
weak. Its ground state is a random singlet (RS) phase 
in which pairs of spins — mostly close together but oc- 
casionally arbitrarily far apart — form singlets. As the 
temperature is lowered, some of these singlets form at 
temperatures of order the typical exchange and become 
inactive. But their neighboring spins will interact weakly 
across them via virtual triplet excitations. At lower tem- 
peratures, such further neighbors can form singlets and 
the process repeats. Concomitantly, the distribution of 
effective coupling strengths broadens rapidly. Eventually, 
singlets form on all length scales and the ground state is 
controlled by an RG fixed point with extremely strong 
disorder: an infinite randomness fixed point. 

This low temperature behavior is in striking contrast 
to that of the pure spin-i AFM chain in which spin-spin 
correlations decay as because of long- wavelength low- 
energy spin- wave (or spinon) modes. In the random sin- 
glet phase, the average correlations decay as a power of 
distance — as — but for a very different reason: A 
typical pair of widely spaced spins will have only exponen- 
tially (in the square root of their separation) small corre- 
lations. But a small fraction, those that form a singlet 
pair, will have correlations of order unity independent of 
their separation; these rare pairs completely dominate the 
average correlations as well as the other low temperature 
properties of the random system. 

Infinite randomness fixed points are ubiquitous in ran- 
dom quantum systems. They probably control phase 
discrete-symmetry breaking transitions in all random 
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FIG. 1. Each connecting line represents a spin-half singlet 
link, (a) Strong randomness spin-| random-singlet phase, (b) 
Low randomness phase of a spin-| chain: valence-bond solid 
-I- spin-^ random singlet. 



quantum systems — in any dimension — g and, in ad- 
dition to the spin-i AFM chain, also control the low- 
temperature properties of a range of random quantum 
phases, jij Because of their ubiquity, further investiga- 
tion of what types of random quantum phases and tran- 
sitions can occur should shed light more generally on the 
combined roles of randomness and quantum fluctuations. 
The simplest cases to analyze are one-dimensional beacuse 
asymptotically exact RG's can be used to extract much of 
the universal low-temperature behavior. In this paper, we 
study random spin-| AFM chains and find that they ex- 
hibit a novel phenomenon: a quantum transition between 
two phases with both phases and the transition governed 
by infinite randomness fixed points. 

We first review what is known about random spin-1 
chains. Pure spin-1 AFM chains behave strikingly differ- 
ently than spin |: their ground state is a non-degenerate 
disordered phase with excitations separated from it by a 
gap. This Haldane gap provides robustness of spin- 
1 chains against weak bounded randomness. [^[0] But for 
strong randomness, spin-1 chains will form a random sin- 
glet phase. As is the case in many random quantum sys- 
tems, there is not a transition directly from the gapped 
phase to the strong randomness phase. Instead, when 
in some local regions the randomness overcomes the gap, 
there will be an intervening region in which there are lo- 
calized gapless excitations but still exponential decay of 
correlations — a Griffiths-McCoy phase. The system un- 
dergoes a quantum transition from this to the random 
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singlet phase as the randomness is increased further. ||Q 
Pure spin- 1 chains with Heisenberg interactions are 
gapless and behave very much hke their spin-i counter- 
parts. 1^,^ We wiU show that random spin-| chains un- 
dergo a phase transition as a function of the randomness 
between two zero-temperature phases: the strong disor- 
der phase is the spin-| analog of the RS phase, with pairs 
of spins forming singlets [Fig. |l|(a)]. Surprisingly, the 
weak randomness phase is also an RS phase, but of an 
effectively spin-i chain superimposed on a Haldane phase 
[Fig. |l](b)]. At a critical disorder, there is a transition 
between these phases, with special behavior at the critical 
point, including a specific combination of spin-i, spin-1, 
and spin- 1 character at low temperatures. 

To make progress, we first review the RG [|l[ analysis of 
random spin-i AFM chains. This proceeds by gradually 
reducing the energy scale, Q. First, the pair of spins with 
the strongest coupling, Jmax — — the initial energy 
scale — forms a singlet [Fig. |2|(a)], and is decimated. 
Virtual triplet excitations cause the two sites neighboring 
the singlet to weakly interact with the effective coupling: 



eff 



JlJr 
Jmax 



(2) 



where Ji, Jr are, respectively, the bonds to the left and 
right of the decimated pair and a = By repeating this 
procedure, we gradually reduce the energy scale, fl, and 
the number of active spins in the chain. In the limit of 
low energy, the random singlet phase emerges and singlets 
form on all length scales. That this occurs for arbitrarily 
weak randomness, as it does, [p^ cannot be convincingly 
shown by this RG as it is initially approximate when the 
distribution of J's is not broad. But its qualitative valid- 
ity for weak randomness is suggested, since Jef / is always 
less than Ji^r due to the prefactor i in Eq. (^. The 
multiplicative structure of Eq. (||) suggests that the dis- 
tribution of J's broadens without bound. This means that 
the perturbative result (^) becomes exact at late stages of 
the RG, 1^ and the universal low energy properties of the 
system can be found exactly. 

The wide distribution of J's allows one to associate 
the renormalized energy scale $7 with the temperature T. 
Bonds stronger than T become frozen, and the remaining 
spins act as though they are free since almost all of their 
couplings are much weaker than T at low temperatures. 

The RG flow is simply parametrized in terms of 



r = In iii 



(3) 



As the RG evolves, ft is reduced, and F increases. At low 
energies the coupling distributions become scale invariant 
functions of /S/F; as F — » oo at the fixed point, the dis- 
tributions become infinitely broad. The density of active 
spins decays as 
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with ■0 = i a universal exponent characterizing the ran- 
dom singlet phase. |2| As V' relates the logarithm of en- 
ergy scales to length scales (1/p), it replaces the exponent 
z which parametrizes power-law energy-length scaling at 
conventional quantum critical points. 

The strong randomness phase of the spin-| chain can 
be understood similarly. Gombining strongly interacting 
neighbors into a singlet yields Eq. (||) with a = |- Strong 
randomness in the J's will guarantee that despite the large 
prefactor (|) the new coupling will almost always obey 
Jeff < Je,r, yielding flow towards the random-singlet 
phase. In Fig. |l|(a), this is indicated by varying length 
triple links representing singlets of spin | . 

When the randomness is weak, the RG for spin | fails 
to reduce the energy scale, suggesting that strong ran- 
domness behavior might not be obtained. To proceed, we 
generalize the method of Monthus, Golinelli and Joliceur. 
j|] Instead of fully decimating strongly coupled pairs of 
spins, we only partially decimate them, eliminating their 
highest energy subspace. Thus, when a spin pair, Sl, Sr, 
is renormalized, its totally ferromagnetic (maximum spin) 
combination is eliminated. This corresponds to breaking 
each spin into spin-i parts — a spin-| consists of three 
spin-i's symmetrized — with each contributing one spin- 
^ to form a spin-half singlet link, leaving a pair of spins 
with S'j^ = Sl, R ~ ^ and modifled couplings between 
them as well as between each one and its other neigh- 
bor. In the ground state, every site must have three links 
joining it to others, e.g., as in Fig. |^. When a link forms 
between two spin-^'s [Fig. ^(a)], both spins disappear and 
the Jf,f f between the remaining neighboring spins is given 
by Eq. (^. As can be seen in Fig. |2|(b), whenever only 
one of an antiferromagnetically coupled pair is spin i , it 
will be decimated, and its partner will form a ferromag- 
netic effective bond across it. Such ferromagnetic (FM) 
bonds can themselves be decimated forming, e.g., a spin 
I from a spin-1 and spin-^ pair; however, no spins greater 
than I can form. We thus see, that as the energy scale is 
lowered, the distribution of effective spins changes. In the 
strong randomness phase, at low energies virtually all the 
active (undecimated) spins have Seff = §■ But this will 
not be the case when the randomness is weak. 

For a spin-1 chain with a narrow distribution of ex- 
changes, i.e. weak randomness, all of the bonds between 
spin-l's would rapidly be partially decimated. The re- 
sulting (approximate) state which has one link connecting 
each site with each of its neighbors is the valence bond 
solid picture of the Haldane phase. [|ll| The scale floo at 
which the last spin is eliminated is the gap. For stronger 
randomness, some double links will form and the gap will 
disappear. But not until a critical randomness is reached 
does the continuous line of links break into finite segments; 
it is this that distinguishes the topological order of the 
Haldane phase from the random singlet phase. 

The phases of a spin-| chain can be understood in a 
related way. With weak randomness, decimation induces 
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FIG. 2. (a) RG rule for the final decimation of two spin 
I's; these are connected by a dashed line which represents a 
spin-i singlet link being formed, (b) Creation of a FM bond 
by the formation of a link marked by a dashed line. Small ar- 
rows indicate preferred relative orientation of the active spins, 
(c) Low energy structure of spin-| chain showing a valence 
bond solid island, composed of effective spin- i's antiferromag- 
netically coupled in its interior, with a spin-1 or a spin-i pair 
ferromagnetically coupled at its ends, separated from other is- 
lands by AFM sea bonds. An AFM islet, made of two spin |'s 
joined by two links, is also shown. Solid arcs represent already 
formed links; the effective spin is noted next to each site. 



singlet links between most neighboring pairs, creating is- 
lands of valence-bond solid. Inside the islands, the active 
degrees of freedom are spin i 's left over from the decima- 
tions with spin I's at the ends of islands [Fig. ||(c)]. The 
islands grow until the entire chain consists of one island 
with only spin i's remaining. At lower energies, these 
spin i's form spin-i random singlets: the ground state 
is thus a spin-i random singlet phase superimposed on a 
(spin-l-like) valence bond sohd; see Fig. |l](b). 

Generally, the low energy structure of a spin-| chain will 
consist of valence bond islands separated by AFM "sea" 
bonds with no links yet formed across them. Each island 
consists of a number — possibly zero — of antiferromag- 
netically coupled active spin-i 's in the interior with each 
end being either spin-1 or two ferromagnetically coupled 
spin-i 's as in Fig. ||(c). The exceptions to this are islets 
consisting of a single AFM bond between two spin-i ends; 
these arise from a pair of spin-| sites connected by two 
links [Fig. |(c)]. There can also be original undecimated 
spin I's. 

It is convenient to describe all this in terms of a purely 
spm-i effective model with a spin 1 represented as a pair 
of spin-i sites with a FM interaction stronger than the 
energy scale, il, and a spin | by an island of three sites 
with two strong FM bonds. This has the advantage that 
coupling distributions and bond types remain indepen- 
dent if they are so initially; thus the number, n, of in- 
ternal spin-i 's in an island is distributed exponentially 
with density oc -B". There are four coupling distribu- 
tions: AFM sea bonds, FM edge bonds, AFM intra-island 
bonds, and (AFM) islet bonds. The other parameters are 



B and q, the fraction of active spins that are in islets. The 
RG flows always broaden without bound the distributions 
of weak (< fl) bonds, [Q justifying the claim that the 
RG is asymptotically exact. In the strong randomness 
limit, i? — > 0, g ^ 0, so that all islands are three spin 
i's strongly FM'ly coupled internally and weakly AFM'ly 
coupled between them equivalent to spin |'s. In contrast, 
for weak randomness at low energies, _B — > 1 and q ^ 0, 
so that an infinite island forms and the system becomes 
equivalent to a random spin-i chain; this then forms a 
spin-i j^s phase. Separating these two zero-temperature 
phases is a novel critical point with non-trivial B and q. 
Both phases and the critical point are controlled by infi- 
nite randomness fixed points. 

To verify the above claims and quantitatively study the 
critical point, we implemented the full RG numerically. 
Initially, H (|l|) is all spin | with the J's uniformly dis- 
tributed in [Jmin, Jmax) and we define 6 = var(ln J). We 
studied 100 realizations of length 5 x 10^, measuring the 
evolution with energy scale of the active spin density p, the 
effective spin distribution, and the coupling distributions. 

For 5 > 5c the chain flows to the Seff = | random- 
singlet phase, while for 5 < 5c flows to the Scff = 
^ random singlet phase. The density p in both random 
singlet phases obeys Eq. (P, as expected, with = ^. 

The critical point is at 5c = 0.22 ± 0.01. The corre- 
sponding fixed point is very different from the stable fixed 
points. The fractions of active spins are (±0.02): 



pi = 0.54 pi = 0.33 PS, ^ 0.13 



(5) 



The appearance of spin-1 excitations may be surprising: in 
pure spin- 1 chains, they do not appear at the ends because 
of the gapless nature of the bulk. At the critical point, the 
active spin density, p, decays with a larger power of F than 
in either phase: 



3.85 ±0.15. 



(6) 
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This implies that the dynamics is faster at the critical 
point than in the adjacent phases. At infinite randomness 
fixed points, '0 also controls the decay of typical correla- 
tions: [Oi 



ln(|(S, -S, 



Cij\i 



J\ 



(7) 



with the random coefficient Cij having a universal distri- 
bution. The average correlations will, however, decay as 
l/|i — at the critical point as in both phases. 

Deviations from the critical point, 6 — 5c, are relevant 
perturbations and grow as F^/"'^- as the energy scale is 
reduced, with the correlation length exponent, we find 



= 1.2 ±0.1 



1^ = 3.2 ±0.3 



(8) 



Many physical quantities are dominated by the almost 
decoupled active spins that remain at scale Q = T cor- 
responding to Tt = hi(fli /T). The magnetization den- 
sity at temperature T and applied field ~ T is the 
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FIG. 3. Evolution of the effective spin fractions as a function 
of r. Squares are spin 1/2, circles spin 1, and triangles spin 
3/2. Filled symbols mark high randomness: 5 — 0.44; empty 
symbols mark low randomness; S = 0.04. The three horizontal 
lines mark the value of the fractions at the fixed point, S = 0.22. 



Ising anisotropy dominates. ||^ But for higher spin, this 
needs exploring. Another intriguing possibihty is three- 
leg spin ladder compounds, if ones can be found with 
combinations of ferro and antiferromagnetic interactions. 
More generally, the model studied here shows how regimes 
with complicated mixtures of effective spins can arise at 
low temperatures from seemingly-simple hamiltonians. 

The spin- 1 AFM chain appears to be the first example 
of a system in which two phases and the transition between 
them are all governed by infinite randomness fixed points. 
How much of this behavior persists in other contexts, in 
particular with lower symmetry or in higher dimensions, 
is a subject for future investigations. 
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Science Foundation via grants DMR-9976621 and DMR- 
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and through Grant No. SFB 484 of the DFG. G.R. would 
like to thank J. P. Sethna for useful discussions. 



the sum of that of the three kinds of spins with weights, 
{Ps{Tt,s)} (see, e.g.. Fig. ||). The hnear susceptibil- 
ity obeys a universal scaling form: x((5, T) w ^ 
T^^^^^M (^{S - 6c)t]/^-''^ /T. For x -> 0, J^{x) ap- 
proaches a nonzero constant, yielding x(T) « l/Tln^^"^" T 
for \S-Sc\ < |lnT|-i/'^-''. For large x, 7V(a;) ~ 
leading to xiS,T) « X{6)/T\n'^T in both random sin- 
glet phases. Near the critical point X{S) vanishes as 
X{d) ~ |5 - J,|(i-2'A<:)i' for \S-6c\ > llnTl-i/V'c^'. Un- 
fortunately, this dip in the susceptibility would be hard 
to observe because of the low temperatures needed. But 
because of the InT in scaling functions, a wide regime of 
the low temperature phase diagram will be governed by 
the critical fixed point with the spin mixtures described 
approximately by the universal fractions in Eq. (^). 

Spatiotemporal correlations can be investigated by neu- 
tron scattering. The magnetic structure factor, S{q,Lu), 
will be dominated at low frequencies by excitations of 
spins that are paired together with energy scale lo. At 
fixed uj, S{q,uj) will show a peak at q ^ Pr((.j)j the typ- 
ical spacing between such spin pairs. At the critical 
point, we also expect some strong ferromagnetic correla- 
tions between widely separated pairs on the same sublat- 
tice. These may give rise to interesting dependence on 6 
of the peak in S{q, to) near the zone boundary. 

The dynamics of nominally pure spin-| Heisenberg 
chains were recently studied experimentally in CsVCl^ 
and CsVBr^, cf. Itoh et al. If mixtures of these, or 
other pairs of compounds can be made with random AFM 
exchange, it should be possible to investigate some of the 
phenomena discussed here. Additional complications that 
would have to be investigated include the effects of random 
anisotropy. For spin-^ random chains, there is consider- 
able robustness of the random singlet-like phases unless 
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